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ABSTRACT 

Planck data towards the galaxy M82 are analyzed in the 70, 100 and 143 GHz bands. A substantial north-south and 
East-West temperature asymmetry is found, extending up to 1° from the galactic center. Being almost frequency- 
independent, these temperature asymmetries are indicative of a Doppler-induced effect regarding the line-of-sight 
dynamics on the halo scale, the ejections from the galactic center and, possibly, even the tidal interaction with M81 
galaxy. The temperature asymmetry thus acts as a model-independent tool to reveal the bulk dynamics in nearby 
edge-on spiral galaxies, like the Sunyaev-Zeldovich effect for clusters of galaxies. 
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1. Introduction 


M82 (NGC 3034, also known as the Cigar galaxy) is the 
largest galaxy of the M81 group in the Ursa Major con¬ 
stellation. It is nearly edge-on (its inclination angle to the 
line of sight is i 80^), and it is the closest galaxy that 
hosts a star burst nucleus, be ing at a distance of only 3.5 
Mpc (Dalcanton et al. [2009). P] Originally classified as a 
dwarf irregular type-II galaxy ^andage fc Brucato||1979 ), 
M82 has been more recently reclassified as a late type SBc 
galaxy after the discovery of a bar (oriented almost along 
the major axis of the galaxy and with about V length) 


and arm structure in its disk (Telesco et al.]|1991 Achter- 
mann & Lacy 1995| Mayya et al.j 2005p. The coordinates 
of the M82 center are R.A. (J2000)= 09^55^52.7®, Dec. 
(J2000)=69^40'46" (corresponding to galactic coordinates 
I = 141.409498^ and b = 40.567579^). It is one of the most 
studi ed galaxies after the Milky Way (see, e.g., Markarian 
1962) and has been extensively observed at all w avelengths. 


fro m the radio band to high energies (see, e.g., 


al. 


Hutton et 


2014p. In the IR, M82 is the brightest galaxy in the 
sEy, a nd shows an outwardly expanding gas ( Hutton et al.j 
2014), which is probably being driven out by the combined 


emerging particle winds of m any stars, together c reating a 
galactic superwind (see also [Gandhi et alT 2011). Stellar 
dynamics measurements of the innermost bulge region of 
the central M82 galaxy show that it harbors a supermas- 
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^ At the M82 distance, 1' corresponds to about 1 kpc. 


sive b lack hole with a mass 3 x 10^ Mq (Gaffney et al. 


1993). 


The M82 rotation curve has been derived in several 
studies using both stellar and g aseous t r acers up t o a ra- 
dius of about 170" 2.9 kpc). [Mayall (1960) and Goetz 


et al. (2014) obtained a rotation curve using optical emis¬ 
sion and absorp tion lines wit h a peak value of about 130 


km s whereas Sofue (1998) used the CO and H I lines to 


find a peak value of about 190 km s“^. Both analyses find a 
steeply rising curve, peaking at about 25" from the galaxy 
center and rapidly declining in a Keplerian fashion beyond 
25", likely indicating that M82’s extended disk mass is miss¬ 
ing, which is v ery peculiar and exc eptional for a disk galaxy. 
More recently, [Greco et al.[ (2012) have presented a K-band 
spectroscopic study and measured the M82 rotation curve 
out to about 4' 4 kpc) on both the eastern and western 

sides of the galaxy. At variance with the previous studies, 
it has been found that the stellar rotation velocity is fiat on 
scales of 1 — 4 kpc (with only a slight decline on the eastern 
side) with a value of ^ 110 — 120 km s“^ out to 4 kpc, thus 
implying an M82 total dynamical mass of ~ 10^^ Mq. 


As shown by De Paolis et al. ([2011 201^ 2015[), CMB 

tne 


data offer a unique opportunity to study the large-scale 
temperature asymmetries far beyond the size typically ac¬ 
cessible with other tools toward nearby astronomical sys¬ 
tems. We study the M82 galaxy by using Planck data 
within a galactocentric radius of 1^. 
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2. Planck data analysis and results 

The sensitivity, angular resolution, and frequency cover¬ 
age of its detectors make Planck a powerful instrument for 
both cosmology and galactic an d extragalactic astrophysics 
( [Planck Collaboration I 2015). Here we use the data of 
Planck 2015 release ([Planck Collaboration I 2015) in the 
bands at 70 GHz of the Low Frequency Instrument (LFI), 
and in the bands at 100 and 143 GHz of the High Frequency 
Instrument (HFI). For a review of Planck results and i nstru¬ 
ment characterist ics, we refer the interested reader to [Buri-[ 
gana et al. ([2013), among others. PlancPs resolution is 13', 
9.6', and 7.V in terms of FHWM at 70, 1 00, and 143 GHz 


bands, respective ly, and frequency maps (Planck Gollabo- 
ration XVI|2015 ) are provided in CMB temperature at res- 
olution corresponding to Nside=‘^04:8 in HEALPix scheme 
(Gorski et al.J|2005[). 


To study the GMB data toward the M82 galaxy, the 
considered region of the sky (shown in Fig. [^in the Planck 
band at 70 GHz) has been divided into four quadrants indi¬ 
cated as Al, A2, A3, and A4. The inner ellipse indicates the 
M82 galaxy extension in the optical band, while the circles 
are at galactocentric radii 15', 30', and 1^, respectively. 

The mean temperature excess (with respect to the mean 
CMB temperature) in jaK in each of the indicated regions 
was obtained in each Planck band at 70^100, and 143 GHz 
with the corresponding standard error [^ As one can see 
from the upper panel of Fig. [^ the A3+A4 region is sys¬ 
tematically hotter than the A1+A2 region with an excess 
temperature almost constant within 15' and 30' in the three 
P/ancA; bands at 70, 100, and 143 GHz. The detected excess 
with respect to the minor axis of the M82 galaxy is therefore 
aligned along the rotational direction, indicating a Doppler 
induced effect modulated by the spin of the galaxy. How¬ 
ever, the temperature asymmetry completely disappears if 
one considers the data within 1^, and this indicates that 
the temperature asymmetry between the two M82 sides is 
substantial only up to about 30 — 40 kpc. 

To test whether the temperature asymmetry we see is 
real or can be explained by a random ffuctuation of the 
CMB signal (which is very patchy, especially on angular 
scales of 0.5^ — 1^), we considered 360 control field regions 
with the same shape as shown in Fig. [^at the same Galactic 
latitude of M82, but at 1^ longitude from each other. For 
each region, we determined the excess temperature profile 
and calculated the average profile and corresponding stan¬ 
dard deviation. The results are shown in the histogram in 
the bottom panel of Fig. As one can see, the temperature 
asymmetry in the 360 control fields in the three considered 
Planck bands is null (within the standard errors), when in¬ 
tegrating within 15', 30', and 1^. Here we note that the 
standard errors within 15' toward the M82 galaxy are quite 
large, as expected, owing to the small number of pixels. 
Therefore, in the three considered bands and within R0.25 
(15'), the temperature asymmetry excess is detected within 
a rather low confidence level (slightly more than 1 standard 
deviation). The situation becomes better when considering 


^ The standard error has been calculated as the standard de¬ 
viation of the excess temperature distribution divided by the 
square root of the pixel number in each region. We have veri¬ 
fied that within the errors, the sigma values calculated in that 
way are consistent with those evaluated by using the covariance 
matrix obtained by a best-fitting procedure with a Gaussian to 
the same distribution. 
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Fig. 1. The Planck field toward the M82 galaxy in the 70GHz 
band. The pixel color gives the temperature excess with respect 
to the mean CMB temperature in jaK. The galactocentric radii 
of the circles are 15', 30' and 1°, respectively. The M82 galaxy 
is indicated by the inner ellipse with major and minor axes of 
10 . 73 ' and 5 . 02 ', respectively. The analysis presented in Section 
2 is performed within the quadrants indicated as Al, A2, A3, 
and A4. The line corresponding to Galactic latitude 40° North 
is also shown. 


the intermediate circle at 30' (R0.50). For example, in this 
case in the 100 GHz band, we detect a temperature excess 
of 50 ± 10 /iK, while the control fields give a temperature 
excess <12 /iK. Therefore, in this region the temperature 
asymmetry is detected with a confidence level of about 3cr. 

Next, we consider the temperature asymmetry with re¬ 
spect to the major axis of the M82 galaxy. In the upper 
panel of Fig. |^we give the CMB temperature asymmetry of 
the A2+A3 regions with respect to the other two regions. 
As in Fig. we also give the temperature asymmetry in 
the 360 control fields, which is consistent with zero in the 
three considered bands and within 15' (R0.25), 30' (R0.50), 
and 1^ (Rl.OO). As one can see, the detected temperature 
asymmetry is substantial (since 40 jjK to 80 /iK) in 
all the Planck bands within at least 1^ (c^ 60 kpc) of the 
M82 center. The very fact that the detected temperature 
asymmetry is almost frequency independent is a clear and 
strong indication of an effect of the galaxy rotation. 

Before closing this section, we stress here that the proce¬ 
dure we followed to test the reliability of the detected tem¬ 
perature asymmetry, that is, of considering the 360 control 
regions, is more reliable than simulating the CMB sky maps 
in each band being considered. While the latter method¬ 
ology of generating sky maps to estimate the error bars is 
mandatory when dealing with the whole sky (as in cosmo¬ 
logical studies), in our case we are considering only fairly 
small regions of the Planck sky maps, and the adopted pro- 
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Fig. 2. Upper panel: the excess temperature toward M82 
in jiK (with the standard errors) of the A3+A4 region with 
respect to the A1+A2 region in the three considered Planck 
bands up to galactocentric radii of 15' (R0.25), 30' (R0.50), and 
1° (Rl.OO). Bottom panel: the same for the 360 control fields 
with the same geometry (shown in Fig. 0 equally spaced at one 
degree distance to each other in Galactic longitude and at the 
same latitude as M82. 


cedure is more reliable since it avoids the ambiguities pos¬ 
sibly involved in the simulations. 

3. Conclusions and further perspectives 

We have considered Planck 2015 release data in the bands 
at 70, 100, and 143 GHz and detected both a east-west 
and a north-south temperature asymmetry toward the M82 
galaxy. We detect a temperature asymmetry of the A3+A4 
versus the A1+A2 regions (that is aligned with respect to 
the M82 rotation axis). This is present in all the consid¬ 
ered Planck bands up to a galactocentric distance of 30 
kpc. Even more robust is the temperature asymmetry of 
the A2+A3 versus the A1+A4 regions, which extends up 
to at least 60 kpc, since it is oriented as the M82 outflow. 

The very fact that the detected temperature asymme¬ 
try is almost frequency independent and is approximately in 
the galaxy spin direction (although probably tilted toward 
the M82 outflow direction) is a rather clear indication of a 
Doppler-induced effect modulated by the spin of the M82 
galaxy. Since the temperature asymmetry is detected far 
beyond the optical extension of the galaxy, it has to bring 
the imprint of the M82 halo, whose rotation axis is proba¬ 
bly tilted in the west direction with respect to the spin axis 
of the disk. This conclusion is reinforced by the observa¬ 
tion that the rotation axis of the ionized emission-line gas 




Fig. 3. Upper panel: the excess temperature in /aK (with 
the standard errors) of the A2+A3 region with respect to the 
A1+A4 region in the three considered Planck bands up to galac¬ 
tocentric radii of 15', 30', and 1° toward the M82 galaxy. Bottom 
panel: the same for the 360 control fields as in Fig. 


is offset from the stellar rotation axis a nd the photomet¬ 
ric majo r axis by an angle of about 12^ (Westmoquette et 


al. 2009). This has been interpreted as possibly resulting 


from the interaction of M82 with the M81 galaxy (which is 
placed approximately in the A3 region in Fig. [^at about 
40' distance from the M82 center) and/or with NGC 3077 
(which lies in the A2 region but beyond the Rl.OO circle) 
in the past. 

We also mention that the rotation of the M82 galaxy 
up to galactocentric distances of a few kpc has also been 
studi ed recently through the kinematics of planetary neb¬ 
ulae ( Johnson et al.||2009 ). The individual radial velocities 
of these planetary nebulae not only agree with the rotation 
curve measured with other methods (see the discussion in 
the Introduction), but moreover, the nebulae at high Galac¬ 
tic latitude (away from the M82 disk) show a clear rotation 
signature, in agreement with our results. 

Before closing, we would like to note that the Planck 
Early Release Compton Source Catalogue has been used, 
together with WMAP (Wilkinso n Microwave Anisotropy 
Probe) data by Peel et al. (2011), to derive the continuum 
spectrum of the M82 galaxy. Then, by a least-squares fit 
to the obtained spectrum by a combination of synchrotron, 
free-free and thermal dust models, the authors have found 
evidence of AME (anomalous microwave emission) from 
this galaxy. However, in a dif ferent way from the analy¬ 
sis presented here, the study by [Peel et al.| (2Qll ) is limited 
to the optical extension of M82 (essentially the 1-2 Planck 

Article number, page 3 of[^ 

















































































pixels of the ellipse at the center of Fig. [B, while our study 
is a “spatial” analysis of the region aroum M82, and, hope¬ 
fully, brings information from its halo, whose presence has 
hitherto been expected. Actually, regardless of the emission 
mechanism that may contribute to the signal at CMB wave¬ 
lengths, the Doppler-induced temperature asymmetry is 
given by 

^ V . . 

— =2reff-sm2 (1) 


where v is the galaxy rotation velocity, i the inclination an¬ 
gle of the galaxy rotation axis with respect to the line of 
sight 1^ and Teff the effective optical depth proper of the 
particular emission mechanism responsible for the temper¬ 
ature asymmetry. For example, for the kSZ case Teff is the 
projected optical depth due to Thomson scattering of th e 
CMB photons on free electrons (see, e.g., Mak et al.|2Qll ), 
while Teff = fS in the case of the cold clouds mechanism 
(where r is the frequency averaged cloud optical depth and 


S is the cloud filling factor, as in |De Paolis et ah 1995). 


As a matter of fact and keeping in mind the object of the 
present study, mechanisms (i), (ii), and (ii;), as detailed 
in footnote 3, require the presence of hot plasma in the 
M82 halo - which is indeed observed up to about 6 kpc 
( Bregman et al.fl995 ) along the minor axis but mig ht also 


be more ex tended as shown by Suzaku observations (Tsuru 
while mechanisms {in) and {v) only involve 


et al 


2007) 


an extended population of cold gas and/or dust, which is 
observed as well 0 

The importance of the methodology we propose, which 
is applied to the M82 case, is that in spite of its simplicity, 
it may allow us to consistently estimate the dynamical mass 
Mdyn (contained within a certain galactocentric distance R) 
of the considered galaxy once the temperature asymmetry 
has been quantified. In fact. 


Mdyn(< R) = ^R = 700 


VTeff smi 


2 

Rkpc Af0 


( 2 ) 


which directly envisages a lower limit to M^yn (obtained for 
Teff = 1) that, in the case of the M82 galaxy, turns out to 
be (AT 80 /iK and T ^ 60 kpc) M^yn > 2.8 x 10^ Mq. 
This is consistent with the observations and can be further 
improved by considering the details regarding the emission 
mechanisms involved (that is estimating Teff). The contri¬ 
bution to the revealed temperature asymmetries can be on 
different scales from the matter outflow from the galactic 
center, the halo rotation, and the stripped matter due to 
M81’s interaction. In general, our method applied to nearby 
edge-on spirals reveals the galactic halo bulk dynamics on 
a rather large scale in a model-independent way. In this 
sense it resembles the SZ effect, which on galactic scales 
cannot work since the e~ temperature is not high enough 


^ These mechanisms are (i) free-free emission, (ii) synchrotron 
emission, (Hi) anomalous microwave emission (AME) from dust 
grains, (iv) the kinetic Sunyaev-Zekdovich (kSZ) effect, and (r') 
cold gas clouds populating the outer gala xy regions (as first pro - 
posed, in the context of the MSI halo, by De Paolis et al.jI995). 
^ In the case of the galactic disks, it can be roughly estimated 
as z = arccos(5/a), where b and a are the minor and major axes 
of the galaxy. 

^ We note that H 2 knots and fila ments in the M82 halo have 
been detected even beyond 5 kpc (|Veilleux et al.||2009 


Beirao 


to produce a substantial effect. The importance of the 
halo’s traced parameters (scale, rotation) by this method 
is also obvious for dark matter and cosmological reasons, 
especially if comple mented by other dynamical and bulk 
motion information (Rauzy & Gurzadyan 1998). Used in 
synergy with other data, this may provide unique keys to 
studying the bulk dynamics, the outflow ejection processes, 
and thus the evolution of galactic systems. 
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